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^ _ We numerically study a one dimensional quasiperiodic system obtained from two dimensional 

■ electrons on the triangular lattice in a uniform magnetic field aided by the multifractal method. 

The phase diagram consists of three phases: two metallic phases and one insulating phase separated 
by critical lines with one bicritical point. Novel transitions between the two metallic phases exist. We 

Q examine the spectra and the wavefunctions along the critical lines. Several types of level statistics 
are obtained. Distributions of the band widths Pb(w) near the origin (in the tail) have a form 
(NJ \ Pb(w) ~ w 13 (P B (w) ~ e~ lw ) (/3,7 > 0), while at the bicritical point P B {w) ~ w~ ' ifi' > 0). 

£S| ■ Also distributions of the level spacings follow an inverse power law P G (s) ~ s~ s (S > ). For 

the wavefunctions at the centers of spectra, scaling exponents and their distribution in terms of the 
Q-/(a)-curve are obtained. The results in the vicinity of critical points are consistent with the phase 
j"^ ■ diagram. 

r/) ' PACS numbers: 71.30.+h, 71.23.Ft, 05.45.Mt 
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I. INTRODUCTION 



A peculiar problem of two-dimensional electrons in a periodic potential with a perpendicular magnetic field has 
' O ■ been attracted much attentions since the Hofstadter butterflji It appears as the spectrum of the underlying one- 
dimensional system called the Harper model 2 which is deduced from electrons on the square lattice in a uniform 
magnetic field . It is also essential in physics of the integer quantum Hall effect^. 

Some of the metal-insulator transitions in one-dimensional quasiperiodic systems have been characterized by mul- 
tifractal structures of band widths and wavefunctions. Seem. A one-dimensional tight-binding model is 



U + iip l+1 + ti-iipi-i + eiipi = Eipi, (1.1) 



where ipi denotes the value of the wavefunction at the i-th site, t% and q are the hopping matrix element and the site 
energy at the i-th site respectively, either or both of them can be taken to be quasiperiodic. The Harper model is the 
t—{ • case where ij = 1 and = A cos(27r ai + 9). When a is an irrational number, it is quasiperiodic. All the eigenstates 
' are extended for A < 2 and are localized for A > 2 with the metal-insulator transitions at A = 2&£. The spectrum has 
a rich structure (the Hofstadter butterfly). The spectrum as well as the eigenstates becomes multifractal 8 . The total 
measure of the bands at the critical point A = 2 is zero with a fractal dimension less than one£. The scaling behaviors 
of the spectra have been extensively studied 8 - by the multifractal analysisi2*ii. Especially the incommensurate limits 
of flux per plaquette, such as the inverse of the golden mean a = have been extensively studied. 

. Recently level statistics of some of the quasiperiodic systems have been investigated and turned out to have the 
behaviors characteristic of criticalitjiiSiii The distributions of the normalized band widths Pb(w)cIw have been 
Q , confirmed that 
O ■ 



P B {w)~<wP (w^O), (1.2) 

and 

P B {w)^ e ^ w (to-foo), (1.3) 

These laws have also been confirmed for a variant of the Harper model at criticalitj*^. A similar type of statistical 
law has been confirmed for the Fibonacci modeiii. Remarkably, the form of the distributions of band widths has a 
similar form as the distributions of the gaps fluctuations observed at the mobility edge of the random systems 15 - 16 . 
Distribution of the energy gaps Pg{s) was also examined. It diverges near the origin and follows an inverse power 

P G (s)~ S - 5 (s^O), (1.4) 



For example, S ~ 1.5 for the critical Harper model(A = 2). 
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One of the aims of this paper is to investigate these quantities for the one dimensional quasiperiodic model obtained 
from two dimensional electrons on the triangular lattice in a uniform magnetic field. This problem was studied by Claro 
and Wannie»i£ , but systematic studies have not been achieved since then. Although the model of two-dimensional 
electrons on the square lattice with next-nearest hopping, which includes the case of triangular lattice as a special 
case, were studied previouslyi*^, statistical techniques such as the multifractal analysis which have been applied for 
other quasiperiodic systems have not been applied to the triangular lattice model. The same model also appears in 
the theory of the junction of three wires of Luttinger liquid^. These situations motivate us to investigate various 
aspects of the quasiperiodic system obtained from the triangular lattice model. 

The organization of this paper is as follows. In Sec|nl we introduce two dimensional electrons on the triangular 
lattice in a uniform magnetic field and obtain the one-dimensional quasiperiodic system. We investigate the Aubry 
and Andre duality^ in this model. In Sec lIIII we investigate the classical orbits of the model and discuss the phase 
diagram. In Sec II VI the distributions of the band widths and the gaps are investigated. In Sec|3 we give a brief 
review of the general formulation of the multifractal analysis and apply it to the spectra and the wavefunctions. We 
confirm the phase diagram conjectured in Sec|n] Sec I VII is the conclusion. 



II. ELECTRONS ON THE TRIANGULAR LATTICE IN A UNIFORM MAGNETIC FIELD 



A. Hamiltonian in real space 



We consider tight- binding electrons on the triangular lattice in a magnetic field (Fig^l. We take the lattice spacing 
to be 1 for simplicity. The Hamiltonian is 

H t a *y C n ^i, m Cn,m 6Xp(^^n+l,m;n,m) t a ^ ^ m C n _|_i jm Q¥L\){iA n 7n . n j r \ rn ^ 

n,rn n,m 

^ ^ C n m ^2Cn,m ^Xp(^^n,m+l;n,m) t b ^ ^ cj t m C n; . m -|_i exp(iv4 n rn ; n5m _)_i) 
n,m n,m 

tc ^ ^ ^ n m ^2Cn+l,m £^p(^-^n,m-l-l;7i-l-l,m) t c ^ ^ C n j^^ ra C n ^ ra J r \ exp(zj4 n ^i_ m ; n rn ^i) 

= H a +H b + H c (III) 

Here t a ,tb and t c are the hopping coefficients for each bond, and c„ jm (c^ m ) is the annihilation (creation) operator 
at site (n,m) : {cjj m ,Ck,i} — ^fc,n<5; m . ^4n,m ; fc,i, fc = n±l, i = m±lisa gauge field on each bond. We impose 
A n ,m;k,i = —Ak.i^n^m so that H to be hermitian. A uniform magnetic field penetrates each triangle with a flux (p = & 
We take the Landau gauge 

Ai+l, m;n,m = 0, A n ,m+l-n,m — 2n(j) and A n , m+ l ;n _|-l !m = 2n<f>{n + -). (H.2) 

thus 2 triangle A ntm ;k,i = §■ A state |*) is written 

|*) = 5>„, m ct im |0). (II.3) 

The Schrodinger equation H\^) = E\ty) is 

- t„(*„-i lTO + ~ h(e 27r ^ n , m _! + e- 2 ^ n # n+1 , m ) 

_f c ( e -2^(»-|)^ n _ 1)m+1 + e 2^(n+^ n+i m) = E y n ^ (n 4) 

We take the form of the wavefunction ^ n . m = e lkym, $ n , then the Schrodinger equation becomes 

- (t„ + t c e- 2 ^ n -^+ lk y )^ n ^ - (t a + t c e 27r ^ ( "+5-^))*„ +1 - 2t b cos(2^n + k y )»£ n = E^ n . (II.5) 

When <fi — ~ (p and q are coprime integers), (III.5|I is periodic with period q. The Bloch theorem tells that one can 
put = exp(ik x n)ip n where i\) n satisfies ip n = 4>n+q, which implies that ^ n +q = e lfex? \I' n . Thus, if we introduce a row 
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FIG. 1: Schematic view of the triangular lattice. 



vector \I/ = ($1, $2, • • • , ^g-i, ^g)' means the transpose of a matrix) and a n (k y ) = t a + i c cxp(27ri^(n + |) — ifcj,) 
and b n (k y ) = 2i{, cos(27r|n + (|II.5|I is reduced to an eigenvalue problem of a finite size matrix 



/ bi(k y ) ai(k y ) 



Hq(k x , ky) 







ai(fc !/ )* b 2 (ky) a 2 (k y ) 

a 2 (fc y )* b3(fc„) ^(fcj,) 

a 3 (k y )* biiky) a^iky) 



e lk ^a (k v )*\ 



a 9 _ 3 (fc y )* b q - 2 {k y ) a q - 2 (k y ) 

a q _ 2 (k y )* b q -i(k y ) a q -i(k y ) 

a 9 _i(fc a )* &<z(fcy) / 



\e- tk *<*a {ky) 
Also, in terms of ip n , ill. 5|) becomes 

- e- ik *(t a + i c e- 2 ^("-^+ lfc ")V„-i - e ik *(t a + t c e 2 "^ n+ ^- lk ^ n+1 - 2i b cos(2^<M + k y )^ n = Eip 
At t c = 0, this is reduced to the Harper equation. We define A = 2|k and fi = ^ then (|II.51 becomes 



1 + ue~ 27 "f (;n-i)+iky 



-1 , 2-rri— (n-\--k —ikv) 

1 + tie i K 2 v ' 



- A COS ( 27T-71 + fc v ) = -E*., 



(II.6) 



(II.7) 



(II.8) 



B. Hamiltonian in momentum space 

We denote k — (k Xy k y ). The electron annihilation operator c(k) in momentum space is 

Cn,m = ,^\ 2 y rffcx y exp^ra + ik y m]c(k). 
The commutation relation for c(k), c^(k) is 

{ C (fe), c^fe')} = (27r) 2 5z(fc x - fc;)^(fc y - fc;), 

where <5z(fc) = J2 n ez ^(^ + 27m). In terms of c(fc), cT(fe), the tight-binding Hamiltonian (|II.1(I is 



(II.9) 



(11.10) 



(11.11) 
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with 

H(k)= ~2t a cosk x c\k)c(k) 

-(t b e- lk « + t c e~ lk * +lk --^)c\k x + k y )c(k x , k y ) 

-{t b e ik v + t c e lk ^ k y- l ^)c\k x - 2tt0, k y )c(k x , k y ). (11.12) 

When cf> = p/q, since k x couples only to k x ± 27T0, we write k x as k x + 2ir<f)j with j £ Z. Here k x is in the magnetic 
Brillouin zone 

- - < k° x < -. (11.13) 
The Hamiltonian -ff (fc) acts on the Hilbert space spanned by 

l*> = E^" ct ( fc ° + 2^,^)10), (11.14) 
=j 

with = ijjj. The Schrodinger equation = E\*f?) is 

- (^e^ + tce-^+^y- 2 ^-^)^-! - (t b e lfc « + i c e^ ,;fc * +2 ™^ + 5))^- +1 - 2t a cos(27i^j + fc°)^ = 

(11.15) 



In (|II.5|I and l|II.15|) , the terms H a and iJb are diagonal respectively. We can also diagonalize H c which is proportional 
to t c by changing the gauge. For example, we take the gauge 

A n+ i, m] , hm = 2ir(f>(n + -), A„ !m+ i ;n , m = 2n(fm, and A„, m+ i ;n+ i, m = 0. (11.16) 

The gauge transformation which transforms from (III.2|I to (|II.16(I is given by 

C'n.m f Cn,m ^Xp(^/n) 



C n,m * C n,m ex P( 



^n,m;n' ,m' ' * -^Ti,m;n / ,m' ~h fn,rn fn',m'i 

f n = 0n(n-l). (11.17) 
In this gauge, the Hamiltonian in momentum space is 

H = 7^2 / dk * I dk y H'(k), (11.18) 

K 2 ™) J-7T 

with 

H'(k) = -2t c cos(k x - k y )c\k)c{k) 

-(t b e- lk y + t a e lk *+^)c) (k x + 2vr0, k y )c(k x , k y ) 

~(t b e lk y + t a e- ik * +ni *)J(k x - 2tt0, k y )c(k x , k y ) (11.19) 

From (|II.19fl . we get the Schrodinger equation : 

- (t b e ik * + t a e- ik *- 2 ™W-^)i> t - X - (t b e- tk y + t a e* k ° +27r ^ e+ ^ e+1 - 2t c cos(27r^ + fc° - k y )^ = E$ t . (11.20) 

Apparently, if one exchange k^ by k x — k y and t a by t c in (|II.15|) . we get (|II.20fl . This is due to the symmetry of the 
triangular lattice. 

C. Duality 

At t c = 0, it is known that l|II.7(l has the duality of Aubry and Andrei who showed the existence of a transition 
between localized and extended states of ipj when (j) is an irrational number. When A > 2, the states are all localized, 
and when A < 2, the states are all extended. At A = 2, all the states are critical. 
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In the present case, take <f> = £ and write 

^=5> 2mM /*- (H.21) 
;=o 

and substitute it into ljll.711 . then 

- (t b e ik » +tc e tk *- tk " +2m ^ l -i'>)f l - 1 - (4e~ 4fe « + tce - *- 4 *"'-^^)/^! -2t a co${2n<t>l + k x )fi = Efi. {11.22) 

When t c = 0, l|11.22|) becomes (I11.15|l by substituting k x — > fc^ an d &j/ — * an< ^ ^ ~~ * X- This is just the Aubry- Andre 
duality when we take the incommensurate limit of cf>. However, when t c ^ 0, i|II.22l) and (|II.15I) are not transformed 
by III. 211) due to the term proportional to t c . 

Because of the symmetry of the triangular lattice, we can consider duality involving t c by putting t a or tb to be 
zero. Let us consider small t a limit. In the limit, ljll.20jl has the Aubry-Andre duality for exchanging tb and t c . This 
implies that there is a duality between A and /i for small t a limit. It relates a state at (A, /i) to the one at (2/i, 4) by 
the transformation l|II.21|) . Thus the phase diagram in (A, (i) should have a localization transitions on the line A = /i 
for small t a limit i.e. large A and /i. 



D. Characteristic Polynomial 



When <f> = p/q is rational, (|II.5p is reduced to the eigenvalue problem of the matrix (|II.6(1 . The eigenvalues are 
determined by the zeroes of the characteristic polynomial 

P(E) =det(E- H q (k x ,k y )). (11.23) 

which has been studied previousljiSiiSiSS. In Refii^, it was shown that P{E) can be written in terms of Chebyshev 
polynomial of order q. In the present case, the characteristic polynomial takes a simple form as follows: 

P(E) = P (E)-Q(k x ,k y ) (11.24) 
Q(k x ,k y ) = {-l)H(tl cos qk x +t q b cos qk y ) + {-lftl COS q(k x -k y ), (11.25) 

where Po(E) is independent of k x and k y . The energy bands are determined by the zeroes of the polynomial P(E) as 
we vary k x and k y . Especially, the edges of energy bands are determined by the minimum and the maximum of the 
function Q(k x , k y ). When k y = 0, they are given by k x = 0, ir. 

In Ref>iS, the total band width W of the triangular lattice is estimated when tb > t a > t c as 

W^(tb-t a )g(q t -^-), (11.26) 



where the scaling function g(x) behaves as 9 ' 3299 when its argument is small. Since t c does not enter the argument, 
the scaling of the total band width of the triangular lattice is the same as the square lattice. This suggests that the 
universality class of the scaling property of the spectral measure of the triangular lattice is the same as that of the 
square lattice. 



III. CLASSICAL ORBITS AND PHASE DIAGRAM 



A. Classical Orbits 



The Hamiltonian (|II.1I) consists of three terms H a , Hb and H c which are noncommutative each other. They are 
diagonalized in different bases as in (|11.5|l . (|11.15|l and (|11.20f) . The "classical" Hamiltonian is thus 

-^classical = 2i a cos k x + 2t b cos k y + 2t c cos(k x - k y ). (III.l) 

In a magnetic field, k y is canonically conjugate to k x and vice versa. Thus to analyze classical orbits, we replace k y 
by x and k x by y. Setting t a = 1, we plot the contour of 

-^classical = COS y + ^ COS X + [I COs(j/ - x) (HI.2) 




^=2.0 M-0 ^3.0^=0.5 




FIG. 3: Contour plots of the classical orbits for (A, fi) = (2.0, 1.0) and (3.0, 0.5). 
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FIG. 4: Phase diagram. In region I and III the wavefunctions (spectra) are extended (absolutely continuous), and in region II 
the wavefunctions are localized (pure points). On the three boundary lines, the wavefunctions (spectra) are critical (singular 
continuous) . 



in Fig. Hfor (X, fi) = (1.0, 0.5), (1.0, 1.0) and (1.0,1.5), and in Fig|3]for (A, fx) = (2.0,1.0) and (3.0,1.5). In FigH 
we see that all the contours for (A, /i) = (1.0,0.5) and (1.0,1.5) are extended in the x-direction and localized in the 
y-direction while, for (A,/x) = (1.0,1.0), there is a separatrix which is extended in both directions. We also see in 
FigC3that the contours for (A,^) = (3.0,0.5) are extended in x-direction and localized in y-direction, while there is 
a separatrix for (X, fi) = (2.0, 1.0). 



B. Phase Diagram for irrational 



From behaviors of the classical orbits shown in the previous section, we may deduce the phase diagram of the 
equation 



- 1 



fie 



-2-Kicp(n— ^)+ik y 



f„_i - 1 + fie 



^2iri(j>(n+ ^ —iky) 



^n+l - X COS (27T071 + k y ) = E^ n . 



(III.3) 



for irrational limit of = 2 . The phase diagram is shown in Fig. |1J One intriguing aspect is the transitions 
between phase I and phase III which are the transitions between metals. Indeed a transition in the quantum case 
is not characterized by an appearance of a separatrix at a certain energy. For the Harper model, it is known that 
metal-insulator transitions occur for whole energies at A = 2. This is generalized to the triangular lattice model we 
consider. 

As an example of an incommensurate limit, in the sections hereafter, we perform numerical scaling analysis for the 
energy spectra and the critical wavefunctions when <f) = ^ approaches the inverse of the golden mean ^ = ^ • A 
standard sequence which corresponds to the continued fraction expansion of — is the Fibonacci series F n , which is 
defined by Fq = Fi = 1, F n = F n -\ + F„_2- F n behaves ~ r™ for large n. By taking p — F n -i,q = F n 

approaches — . F n is called a Fibonacci number and n is referred to Fibonacci index. 

To take this incommensurate limit of l|IH.3|) . the off-diagonal terms in ljIII.3|) need a special care. Namely, when 
|U = 1, these terms can be zero if the exponential becomes —1. For the sequence above, this actually happens when 
q = F n with n = 3£ + 1 for some integer £. In that case, the energy spectrum has no dependence on k x , and the 
dispersion relation is flat. 



— £ 
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IV. LEVEL STATISTICS 



Consider 



when k y = 0. On the critical lines, the spectral measure and the wavefunctions are expected to 



show characteristic behaviors of criticality. See Figj4] In order to obtain the distributions of the band widths, the 
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FIG. 5: Distributions of the band widths for (A, fi) = (2.0, 0.4). 



gxij matrices (|II.6(1 are diagonalized. The normalizations are 

poo 

/ P B {w)dw = 1 
Jo 



(w) 



wPb{w)cIw 



Similarly the distributions of the gaps Pg(s) are obtained and normalized by 



poo 

/ P G {s)ds = 1 
Jo 

poo 

/ sP G (s)ds = 1. 
Jo 



(IV.l) 



(IV.2) 



As we discussed in Sec Mil the edges of the energy bands are found at k x — 0, n and k v — 0, n. Thus, to study the 
measure of the spectrum of (|III.3() . it is sufficient to study those points in the Brillouin zone. When (j> = £ is a rational 



number, the problem is reduced to the eigenvalue problem of the finite size matrix (|II.6() . Furthermore, when q is 
odd and k x = 0, n and k y = 0, tt, the matrix - 6|) reduces to a tridiagonal form by the symmetric and antisymmetric 
eigenstateA 

For [l = 1 with q = F n and n = 3£ + 1 , as we noted above, the hopping term becomes zero at a bond and all the 
band has zero width. Thus we study only the case of q = F n with n = 3£ when fj, = 1. 



A. Distributions of Band Widths 



Consider the distributions of the band widths along the line A = 2 with /i = 0.2, 0.4, 0.6 and 0.8. In Fig. \S\ Pb(w) 
at (A, /i) = (2.0, 0.4) for q = F n with n = 25, 27 and 28 are plotted. It shows convergence to a limit, indicating the 
existence of a limit of the distributions of the gaps for the incommensurate flux ip. For < [i < 1 the distributions 
depend on [i. The semi-log plots of Pb{w) is shown in Fig. |H| One sees the linear behaviors for large w, implying an 
asymptotic form 

P B (w)~ e -'i w , as w-^oo, (IV.3) 

where 7 > 0. The optimized values of 7 are shown in Table [I| for several /i's. The inset of FigJS] shows Pb{w) near 
the origin which indicates that the distributions of the band widths Pg(w) are zero at the origin with a power law 
decay. To characterize this behavior, we make an ansatz 



Pb(w) ~ w , as w — * 0. 



(IV.4) 
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0.5 1 1.5 2 2.5 3 3.5 

W 

FIG. 6: Semi-log plots of the distributions of the band widths for (A,/i) = (2.0, 0.4). Inset: Log-log plots of the distributions 
of the band widths for (A, /i) = (2.0,0.4). As seen in the inset, the convergence of the distributions near zero at large n is 
relatively slow. For each value of (A,/x), we choose a stable part of the distributions to obtain the exponents. This procedure 
potentially underestimates the values. 




0.5 1 1.5 2 2.5 3 3.5 

W 

FIG. 7: Semi-log plots of the distributions of the band widths for (A,/i) = (1.0, 1.0). Inset: Log-log plots of the distributions 
of the band widths for (A,/x) = (1.0, 1.0). 



where (3 > 0. The optimized values of (3 are shown in Table [i] One sees that f3 becomes smaller as approaching to 
M = l. 

Next, we investigate Pb(w) on the other lines /x = 1 and A = 2fi. In Fig. [7| , the semi-log and the log- log plots of 
Pb(w) are shown for (A,/x) = (1.0, 1.0). We find similar type of behaviors l)IV.3jl and 1IV.4|I for the A = 2 line. We 
also investigate the critical line A = 2/i and find similar type of behaviors. We collect the values of /3,7 in Table [I] 

These behaviors of Pb(w) on these lines are consistent with the behavior of Pb(w) in other quasiperiodic 
mode l 12 ^?' 14 thus gives a support for the phase diagram of FigEl 

B. Distributions of Gaps 

The distribution of the gaps at (A, fi) = (2.0, 0.0) has been known to follow an inverse power lawii which diverges 
at the origin 

P G (s) ~ s- 5 (IV.5) 

with 5 - 1.5. In Fig. |U the distributions of the gaps for (A, (i) = (2.0, 0.4), (1.0, 1.0) and (2.0, 1.0) are shown. It is 
clear that Pg(s) shows a power law of the inverse. The estimated value of S is ~ 1.5 for these cases. We also investigate 
other points on the lines A = 2, fj, = 1 and X — 2fi and find a similar behavior with <5 ~ 1.5 within statistical error. 
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FIG. 8: Log-log plots of the distributions of the gaps for (A, ft) = (2.0, 0.4), (1.0, 1.0). 



This behavior of Pg(s) shows that the spectra are singular continuous on these lines, gives a further support for the 
phase diagram FigJI] Also the value 6 ~ 1.5 seems to be a characteristic quantity for this model. 



C. Bicritical Point 



We also investigate the point (A, /x) = (2.0, 1.0). The result is shown in FigJ5| In sharp contrast to other points on 
the critical lines, it shows the inverse power law 

PbH-u)-"', (IV.6) 

(/?' > 0) for whole the range. We estimate the exponent of the law as (3' ~ 1.4. This implies that the spectrum at 
this point is a qualitively different fractal-like set. On the other hand, the gap distribution is shown in Fig llOl It is 
an inverse power law 

P G (s) ~ s- s , (IV.7) 

(5 > 0) with the exponent S ~ 1.5, analogous to the ones found for other critical points. Thus the band width 
distribution gives a finer characterization of the spectra than the gap distribution. 



V. MULTIFRACTAL ANALYSIS 



We apply the method of multifractal analysis^ to the spectra and the critical wavefunctions. In Ref the entropy 
function was introduced which reformulates the theory along the way that standard statistical mechanics is formulated. 
We use it in our analysis. 
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FIG. 9: Log-log plots of the distributions of the band widths for (A,/i) = (2.0, 1.0). 
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FIG. 10: Log-log plots of the distributions of the gaps for (A, fj) = (2.0, 1.0). 

A. Review of Multifractal Analysis 

We consider quantities Zj and their probability measure of a fractal- like set. Though we shall only consider 
the cases where U or pi is a constant, a general formulation is reviewed for convenience. It is natural to consider 



A fi /3 7 8 

2.0 0.0 2.5 1.4 1.5 

2.0 0.2 2.5 1.3 1.5 

2.0 0.4 2.5 1.2 1.5 

2.0 0.6 2.3 1.2 1.5 

2.0 0.8 2.1 0.9 1.5 

2.5 1.25 2.1 0.9 1.5 

3.0 1.5 2.3 1.1 1.5 

4.0 2.0 2.4 1.2 1.5 

0.0 1.0 2.6 1.6 1.5 

0.5 1.0 2.6 1.5 1.5 

1.0 1.0 2.4 1.3 1.5 

1.5 1.0 2.2 1.0 1.5 



TABLE I: Estimated exponents on the critical lines. For definitions of /3, 7 and <5, see Eqs. lIV.3ll . Ill V .41 . and llV.5t respectively. 
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distributions of logarithm of h 

Inl 



-, i.e. U = exp(— ne). (V.l) 
n 

As n becomes large, li approaches zero , but ej takes a hnite nonzero value for critical points. We introduce a scale 
index a, as the exponent of measured by U as 

Pi = if*, ai = ---ln Pi . (V.2) 
£j n 

We write the number of li whose scale index lies between e and e + de, a and a + da as fi(e:, ajdeda. We take an 
ansatz that 17 (e, a) has the following scaling form for large n 

Q(s,a) = exp[nQ(e,a)], (V.3) 



where Q(s, a) can be seen as a kind of entropy function. 

Followingifliii, we consider the generalized partition function 



The generalized free energy is 



Using Q(e, a), IV. 5|) is written 



T(q,0) = £>?Zf (V.4) 

i 

= exp[-ngi(aig + /?)]. (V.5) 



G(g,/3) = ilnr( g ,/3), (V.6) 



T(q,/3) = J de J daexp[n[Q(e,a) - (aq + /3)e]]. (V.7) 



For large rt, the maximum of the exponent dominates the integral and gives 

G(q, /?) = Q«£>, (a)) - «a)<z + /3)(e), (V.8) 
where (e) and (a) give the maximum of Q(e,a) — (aq + [3)e, so we have 

dQ(e,a) 



de \e=(e), a =(a) 
and 

dQ(e,a) 
da 



(a) + 0. (V.9) 



e)q. (V.10) 



Thus G(q, (3) is obtained from Q(e, a) using (|V.8()I|V.9|) and (|V.10|1 . From 1V.10|) . the maximum of Q(e,a) with 
respect to a occurs when q = 0. On the other hand, once G(q,f3) is calculated, (e), (a) and <5((e), (a)) are given by 

(e) = ~G{q,p), (a)(e) = ~G(q,/3), (V.ll) 

and 

<2((£>>W) =G{q,P) ~Q — ^ P — gp — ■ (V.12) 

Since (e) and (a) are functions of q and f3, different regions with scaling indices e and a are explored by changing the 
values of the parameters q and (3. Thus Q({s), (a)) is implicitly a function of q and (3. 
The limit of G(q, (3) for large n, may be obtained by 

G(q,0 B (q))=O, (V.13) 



13 



and P c (q) can be regarded as a set of generalized dimensions. The scaling index (e) c which corresponds to /3 c (q) could 
be considered as being a representative for a particular value of q. 

From fVJ^l JvTflll and fvT^ . we see that Q{(s), (a)) at the critical point satisfies the relation 



where f((a) c ) is given by 



Q«e>c»c) = ^|^U< e ) c <£} c = /(W<£> c , (V.14) 



/((«) c ) = q- £ l B= ( e ) e - (V.15) 



By substituting (|V.15f> into ljV.9|l and (|V.10fl . we obtain 

f({a) c ) = (a) c q + f3 c (q) (V.16) 



and 



respectively. And (|V.16fl and (|V.17(I give 



d f(( a )c) i rvl7 ^ 
9 ' (V ' 17) 



<«> e = (V.18) 

Thus once (3 c (q) is known by solving (|V.13(I . (a) c and f{{ a) r ) are obtained from fVU6|l and (fVTTSji . In terms of /(a), 
the density function of £ and a f2(e,a) is written, using (|V.3|I and (|V.14|I as 

Q((e), (a) c ) = exp[n( £ ) c /({a) c )] = (f>J /( < a >.\ (V.19) 

where (Z) c = exp(— n(e) c ) is a representative length, and (e) c and (a) c are functions of q [see l|V.ll|) and l|V.18|) ]. 
f(a) can be considered to be a set of generalized dimensions. In numerical approach, we calculate f(a) for a given 
fractal-like object for a finite n and extrapolate it to the limit n — > oo. From G(g, j3) l|V.6|) . G n (q, (3) at finite n should 
behave as G n (q,f3) ~ G(q,f3) + O(-). Thus we should extrapolate G n (q,/3) as a function of i and estimate the limit 
for n — > oo. 

We denote the support of /(a) by [a m i n ,a max ] and the value of a which gives the maximum of f(a) by «o- 

B. Spectrum 

In this section, scaling properties of energy spectra are analyzed by the multifractal analysis. 

1. Multifractal Analysis of Spectrum 

We apply the general formulation to characterization of the energy spectra. Take the band widths as variables 
and let pi be constants 

Pi = y ~ ^ ( v - 2 °) 



where r is the golden mean. As seen from (|V.6|) . 

G(q,f3) = -qln T + F((3), = G(0,f3) = -lnr(0,/3). (V.21) 

n 

For each q, the critical value (3 is determined and vice versa. From IjV.llj) and (|V.21fl . the scaling index (e) is witten 



14 



0.5 



0.4 



0.3 
0.2 



0.1 





0.4 0.42 0.44 0.46 0.48 0.5 0.52 0.54 0.56 

a 

FIG. 11: a-f(a) curves of the spectrum for = (1.0, 1.0). 

Thus (e) depends only on f3. From l|V.ll|) and i|V.21|) . (a) is related to (e) by 

(a)=JnT/(e), (V.23) 

and a) is nonzero only for a satisfying l|V.23|l . thus depends only on e. We may put Q(e, a) as S(e), and from 
EH} and (fVT4)l we get 

S(e)=F(f3)+eP = ef(a). (V.24) 

Then /(a) is calculated from the formula 

f(a) = ^-. (V.25) 
2. Numerical Results 

We apply the method above to the spectra of our model. In Fig JIll we show a-f (a) curve for the spectrum 
at (A, /Lt) = (1.0,1.0). The estimated values of a m { n , a max and ao are 0.421, 0.547 and 0.495 respectively. These 
values coincide with the corresponding values of the Harper model^. We also investigated other points on lines 
A = 2, fi = 1, A = 2/i in Figj4j Except for (A, (i) = (2.0, 1.0), it turns out that the estimated values of a m i n , a max and 
ao are the same as those of the Harper model. This implies that the universality class for these lines is the same as 
the Harper model. This is consistent with the scaling of the total band widths l|II.26(l for A = 2, fj, < 1. On the other- 
hand, a-f (a) curve for the spectrum at (A,/i) = (2.0,1.0) has a different shape as shown in Fig J12l We see that 
c*min = 0.381, a max = 0.755 and «o = 0.498. Thus the universality of this point is different from the Harper model. 

C. Wavefunctions 

We investigate scaling properties of the wavefunctions by the multifractal analysis. We concentrate on the eigen- 
functions at the centers of the spectra. This enables us to confirm the phase diagram Fig0] 

1. Multifractal Analysis of Wavefunctions 

We apply the general formulation to characterize the wavefunctions. Take squares modulus of the wavefunctions 
to be variables pi, while take U to be constants 
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FIG. 12: a-f(a) curves of the spectrum for (A,/i) = (2.0, 1.0). 



Thus e is a constant in this case. From l|V.5(l and (|V.6J| . one has 

G(q,0) = -/3e + G(q,O). 
Using IjV.lip and i|V.12|) . we obtain the generalized entropy 

Q(e,a) = G(q,0)-qe(a), (a) = -\^§^-- 
Since Q(e, a) is nonzero only for (a), we write it as S'(a). From l|V.14|) . we have 

/(«) = 



(V.27) 



(V.28) 



(V.29) 



We calculate f{a) for finite Fibonacci index n by this formula and extrapolate them to n — > oo. Actually only a 
part of f(a) is required to distinguish localized, extended and critical states. For a localized state, /(a) has a point 
support and takes nonzero value only at a mm = and a max = oo and f(a min ) = 0, /(a m ax) = 1- For an extended 
state, it has a m j n = OL max = 1 and f(a m - ln ) — /(a max ) = 1. For a critical state, f(a) may have a finite interval 
[a m in, a m ax] as a support and /(a) takes various values. We shall use this method to distinguish states near critical 
points. 



2. Numerical Results 



We numerically obtain the wavefunctions at the centers of the spectra for odd q on the A = 2, fi = 1 and A = 2/j, 
lines in the phase diagram Figj4] For [i = 1, we investigate q = F n with n = 3£ as well as 3^+1. Although the 
dispersion relations are flat when n, = 3£ + 1, we find that the wavefunctions still show a characteristic behavior of a 
critical state. 

In Fig J13l the square moduli of the wavefunctions at the centers of spectra for (A,/i) = (1.0, 1.1), (1.0, 1.0) and 
(1.0,0.9) are displayed for n = 21 and F n — 17711. From these figures, we see that the wavefunctions are extended 
for (1.0, 1.1) and for (1.0,0.9), and critical for (1.0, 1.0) which is in accord with the phase diagram FigQ] In Fig J14l 
the square modulus of the wavefunctions at the band center for (A,/i) = (2.0, 1.0), (2.0, 1.1), (1.9, 0.9) and (2.1,1.0) 
i.e. in the vicinity of the bicritical point of FigJ3]are displayed for n = 21 and F n = 17711 (n = 18 and F n = 4181 
for (2.1, 1.0)). It is rather clear that the wavefunction is extended for (2.0, 1.1) and (1.9,0.9), localized for (1.0,0.9) 
and critical for (2.0, 1.0). To draw convincing conclusions, however, it is necessary to study the scaling properties by 
multifractal analysis. 

We plot a min for = (1.0, 1.1), (1.0, 1.0) and (1.0,0.9) in FigHU For (A, fi) = (1.0,0.9) and (1.0,1.1), it 

is clearly seen that a m ; n extrapolates to 1 for n — > oo. On the other hand, a m i n is extrapolated to 0.358 for 
(A,/z) = (1.0, 1.0). This value of a m ; n is actually the same as the one found in the Harper model 5 within statistical 
error. As shown in FiglTHl f(a m j a ) extrapolates to 1 for (A, fi) = (1.0,1.1) and (1.0,0.9), and for (1.0,1.0). The 
behaviors of a m i n and /(a m i n ) in Fig J15N16l indicate that the state is extended for (A,/i) = (1.0, 1.1) and (1.0,0.9), 
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and critical for (1.0, 1.0). This confirms a part of the phase diagram in FigJI] especially the metal-metal transitions 
at \i = 1.0. 

In Fig. 1171 a m in's are shown for the states near the bicritical point (A, /i) = (2.0, 1.0), (2.0, 1.1), (1.9, 0.9) and 
(2.1, 1.0). Also Fig. H8l shows /(a m i n )'s for them. For (2.0, 1.1) and (1.9, 0.9), both a min and /(a m i n ) are extrapolated 
to 1, telling that the state is extended. For (2.1, 1.0), both a m i n and /(a m in) are extrapolated to 0, which means that 
the state is localized. At (2.0, 1.0), the convergence of a m i n seems slow but the plots show a tendency to converge to 
a finite value near 0.47. Similarly f(a m [ n ) converges to zero. Thus we conclude that the states are critical at (2.0, 1.0) 
and it is the bicritical point of metal-insulator and metal-metal transitions. See FigJI] 

Next Fig J19l and Fig J20l show the scaling of a m [ n and /(a m i n ) respectively near (A, fi) = (3.0, 1.5). The extrapolated 
values of a min and /(a m ; n ) are consistent with the phases diagram FigJI] We investigate a m i n and /(a m ; n ) for other 
points on the critical lines and the results are consistent with the (A, /i)- phase diagram FigJI] 

Let us next turn to whole a-f(a) curve. In Fig l21l the a-f(a) curve at (A, /.i) = (2.0, 0.4) is shown. Within statistical 
error, a m i n is 0.358. This value of a m i n holds on the critical lines except in the vicinity of the bicritical point, where 
a m in slightly changes about 0.05. In Fig J21l the value of ao which gives the maximum of /(a) is observed to be 1.31 
which is also the same as the Harper models. It is the same for the other points on the critical lines A = 2.0 and 
A = 2/i within statistical error, except in the vicinity of the bicritical point where ao slightly changes within 0.05. We 
may interpret these slight changes of a m i n and ao in the vicinity of the bicritical point as an effect of slow convergence 
there (see below). For /(a) for a > ao, the convergence of /(a) is not good, especially for a max . 

For the critical line fi — 1, ao splits for n = 2>£ and n — 3^+ 1. In Fig l22l plots of ao and a m i n for (A, /i) = (1.0, 1.0) 
are shown. It is seen that ao goes to different values for n — Zt and n = 2>l + 1. Similar behavior of ao is observed 
for other points on the fj, = 1.0 critical line. This implies that the universality class is different for these two series. 
Also the a-f (a) curves for n = 3£ and 3£ + 1 are shown in Fig J23l The estimated values of a mm and ao for n = 3£ 
and 3£ + 1 with /x = 1 are shown in Table [H] . 

For the bicritical point (A, fi) — (2.0, 1.0), the split of 3£ and 3£ + 1 is not so obvious in the numerical data. For 
example, Fig. 1171 shows a m i n for n = 12, 13, • • • and they don't clearly split into two series. Also the convergence is 
not as good as those for other points. The numeric for (A, /i) — (2.0, 1.0) in Table ITT1 is obtained based on both series. 
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FIG. 14: Plots of the wavefunctions at the centers of the spectra for (a) (A, fx) = (2.0, 1.0) , (b) (2.0, 1.1), (c) (1.9,0.9) and (d) 
(2.1, 1.0). Here q = 17711 for (a), (b) and (c), and q = 4181 for (d). 
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FIG. 15: Plots of amir, vs. i near (A, fi) = (1.0, 1.0) with n = 12, 13, 15, 16, 17, 18, 19, 21 and 22. 



VI. CONCLUSIONS 

We study two dimensional electrons on the triangular lattice in a uniform magnetic field and the one dimensional 
quasiperiodic system obtained from it. We conjectured a phase diagram of the one dimensional model as in Fig. |4] 
As a typical example, we investigated the incommensurate limit of the golden mean via the level statistics, namely the 
distributions of the band widths and the gaps, and scaling properties of spectra and wavefunctions on the conjectured 
critical lines. For level statistics, we find the characteristic behaviors similar to the ones previously found for other 
quasiperiodic models. We also obtain a-f (a) curve for spectra and the wavefunctions at the centers of the spectra. 
As for the spectra, a-f (a) curve is the same as one in the Harper model on the critical lines except for the bicritical 
point. For the wavefunctions, we find that a m i n is the same as the Harper model except near the bicritical point. For 
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FIG. 16: Plots of /(Qmin) against ^ near (A,/n) = (1.0, 1.0). 
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FIG. 17: Plots of Q min against ± near (A, /i) = (2.0, 1.0). 



the line A = 2, a-f(a) curves are the same as the Harper model. For fj,=l, the Fibonacci sequence F n splits into two 
classes n — 31 and n — 31 + 1 according to the appearance of zero of the hopping terms. The dispersion relation is 
flat for n = 31 + 1. Also the a-f (a) curve of the wavefunction is different for n — 3l and n = 3£ + 1. At the bicritical 
point where the triangular lattice symmetry retains, both level statistics and multifractal analysis show qualitivcly 
different behaviors from those of other critical points. 
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FIG. 18: Plots of /(Q min ) vs. i near (X,n) = (2.0, 1.0). 
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FIG. 20: Plots of /(Qmin) against ± near (A,//) = (3.0, 1.5). 
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